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ABSTRACT 

Cell mechanical functions like locomotion, contraction and 
division are controlled by the cytoskeleton, a dynamic 
biopolymer network whose mechanical properties remain 
poorly understood. We perform single-cell uniaxial stretch- 
ing experiments on 3T3 fibroblasts. By superimposing 
small amplitude oscillations on a mechanically prestressed 
cell, we find a transition from linear viscoelastic behavior 
to power-law stress stiffening. Data from different cells 
over several stress decades can be uniquely scaled to ob- 
tain a master-relation between the viscoelastic moduli and 
the average force. Remarkably, this relation holds indepen- 
dently of deformation history, adhesion biochemistry, and 
intensity of active contraction. In particular, it is irrele- 
vant whether force is actively generated by the cell or exter- 
nally imposed by stretching. We propose that the master- 
relation reflects the mechanical behavior of the force bear- 
ing actin cytoskeleton, in agreement with stress stiffening 
known from semiflexible filament networks. 

INTRODUCTION 

Mechanical forces are essential for biological systems and their 
interaction with the environment. Tissues can generate forces, 
but forces also influence tissue development, as in embryogen- 
esis, bone growth or scar tissue formation In vitro, single 
eucaryotic cells develop tension spontaneously on a substrate 
a, and respond to extracellular elasticity ft. Cell mechan- 



ics is controlled by the cytoskeleton a highly dynamic 

biopolymer network far from thermodynamic equilibrium. Its 
main components are actin, tubulin and intermediate filaments. 
Different types of crosslinking proteins join these filaments into 
a myriad of geometries. Particularly interesting for cell me- 
chanical properties are actin filaments, which can generate ac- 
tive forces through interaction with myosin motor proteins. 

In the last years, much progress has been made in the study 
of reconstituted biopolymer systems with a known and lim- 
ited number of components. The rheological behavior of such 
biopolymer networks has been experimentally studied and ex- 
plained by simple models (0, lioh - Relating the so 
gained framework to living cells, however, is a daunting task. 
Unlike passive networks, the actin cytoskeleton actively gener- 
ates forces. Further, because of cytoskeleton restructuring or 
changes in motor activity, the mechanical behavior evolves in 
time. This evolution is controlled by signaling cascades, which 
in turn are influenced by mechanical perturbations. 

Probing the mechanical properties of the actin network at the 
micrometer scale has become possible thanks to newly devel- 
oped microrheological techniques lillUl2l.ll3l.ll4l) . However, 
cell heterogeneity complicates the extrapolation of local prop- 
erties to the whole cell scale. Many-cell experiments, such as 
those on tissue or cell-populated gels (Il5lll6lll7l) . are efficient 
and reproducible, but they average over heterogeneous sets of 
cells interconnected by extracellular matrix. This hampers an 
interpretation in single cell terms. The quantitative study of the 
mechanical properties of entire, single cells is thus an essential 
step towards a global understanding of cell mechanics. 

We report uniaxial stretching experiments on single 3T3 fi- 
broblasts, suspended between two parallel, biochemically func- 
tionalized walls (see Fig.Q. This microplate rheometer allows 
us to clearly distinguish active behavior from the passive me- 
chanical response to applied forces or strains ( 18). 
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TOP VIEW 




Figure 1: Image of an actively contracting fibroblast between two 
fibronectin coated microplates. The distance between plates is ~ 5fim. 



MATERIALS AND METHODS 
Experimental setup 

We describe the main features of our cell-rheometer, schemati- 
cally shown in Fig.[2] This is an improved version of the orig- 
inal, home built micromanipulation setup previously described 
3T3 Fibroblasts (19, 20) are held between two parallel, 
coated glass microplates. One of them is rigid and provides a 
reference point. The other microplate is thin and narrow and 
therefore flexible, with an effective stiffness ~ 100 nN//im, 
which is of the order of the elastic modulus of a typical fibrob- 
last. Its bending provides the force acting on the cell in the y- 
axis (see Fig.|5J. The flexible microplate is translated by means 
of a piezoelectric actuator P-841.40 (Physik Instrumente, Karl- 
sruhe, Germany) with a resolution of 1 nm. An optical fiber in 
contact with the flexible microplate couples microplate bending 
to translation of the emergent laser beam. The optical fiber is 
etched with hydrofluoric acid to a diameter of 6 /im, so that its 
stiffness is far lower than that of the flexible plate and does not 
interfere with the force measurement. The cell is illuminated 
with green light and observed with an Axiovert 135 microscope 
(Zeiss, Oberkochen, Germany). A dichroic mirror separates 
the green illumination light from the He-Ne laser beam, which 
reaches a S-1880 two dimensional position sensitive detector 
(Hamamatsu photonics, Japan) through one of the microscope 
ports. A personal computer reads the signal from the position 
detector, controls the piezoelectric actuator and calculates the 
normal force F and the cell length I. The precision in the mea- 
surement of the fiber position is about 100 nm. In most cases a 
high resolution in time is not needed and averaging can improve 
the resolution to 30 nm. By controlling the piezoelectric trans- 
lator, a feedback loop can impose user defined force or length 
histories. The response time of the piezoelectric translator lim- 
its dynamic measurements to frequencies below 30 Hz. Due to 
the large length of the flexible microplate, its tip deflects by less 
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Figure 2: Schematic of the micromanipulation set-up. A fibroblast is 
held between two coated microplates. The deformation of the flexible 
microplate gives the force F acting on the cell. The position of the 
laser beam emerging from the optical fiber which is in contact with the 
tip of the flexible plate is detected using a position sensitive detector 
(PSD). A personal computer reads the signal from the detector and 
controls the piezoelectric translator. 



than 6 arc min during an experiment. Hence, the experimental 
geometry can be described as two parallel walls, which can be 
separated by a translation in the perpendicular direction. 

Temperature Control. The temperature of the cell-culture 
medium inside the chamber is controlled by means of two ITO- 
coated glass slides, one on the bottom of the chamber, the other 
one above, with holes for the microplates and the optical fiber 
to pass through. In order to avoid convection in the medium, 
these slides are kept at different temperatures, imposing a tem- 
perature gradient pointing upwards with a magnitude of ~ 1 
°C/cm. 

Microplates. The glass microplates used for the experi- 
ment are obtained by pulling glass strips (Vitrocom, NJ) as de- 
scribed previously (18), using a modified P-97 Flaming/Brown 
micropipette puller (Sutter Instruments, Novato, CA). Mi- 
croplates are then cleaned and coated with fibronectin from 
bovine plasma (Sigma-Aldrich) or 3-aminopropyl triethoxysi- 
lane (Sigma-Aldrich) and glutaraldehyde (Fluka Chemie, 
Buchs, Switzerland). 

Cell culture. Experiments where the cells stick to 
fibronectin-coated plates are performed in ISCOVE medium, 
with 25 mM HEPES buffer and 10% fetal bovine serum (FBS). 
Experiments with a glutaraldehyde coating start in pure saline 
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solution (HBSS) to avoid inactivation of the coating by proteins 
and amino-acids. After attachment of the cell to both plates, 
FBS is added to a final concentration of 2%. All cell culture 
reagents are from Gibco (Invitrogen, Carlsbad, CA). 3T3 fi- 
broblasts are obtained from the German Collection of Microor- 
ganisms and Cell Cultures (DSMZ, Braunschweig, Germany) 



RESULTS 

Active response at constant cell-length 

We first perform experiments to characterize the response of 
fibroblasts to their presence in the rheometer (see Fig.0. 
To stimulate contractility, we use high serum concentrations 
of 10%, and fibronectin mediated adhesion using coated mi- 
croplates. Fibronectin binding to integrins is known to trigger 
the formation of focal complexes, which connect the extracel- 
lular matrix to the actin cytoskeleton (I22ll23h . To minimize the 
mechanical perturbation to the cell, we keep the cell length £, 
given by the distance between the rheometer microplates con- 
stant and measure the force F. A typical response is shown in 
Fig. |3] After contact with the fibronectin coatings, most cells 
generate contractile tension, the force F reaches values up to 
1 /iN and eventually decays to zero. The force relaxation is an 
adaptation to the constant length condition, since active con- 
traction can be induced again by a sudden change in cell length 
I. The behavior is reproducible only in its broad, qualitative 
features. The force and time scales are strongly cell-dependent. 
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Figure 3: Force as a function of time at constant cell extension, 
recorded immediately after incorporation of the fibroblast into the 
rheometer. No significant cell shape alterations are seen throughout. 
The experiment is performed using fibronectin mediated adhesion. 



Dynamic viscoelastic moduli 

In order to probe cell rheological properties during complex 
adaptive processes, we superimpose sinusoidal length oscilla- 
tions at small amplitudes and high frequencies. These are cho- 
sen so that the corresponding maximum rate of change of force 
is at least two orders of magnitude above the values typically 
seen in active contraction. 

Cells have both an elastic and a viscous response to imposed 
length oscillations qjj, l24f) . As the viscous response depends 
on deformation rate instead of absolute deformation, it causes 
a phase-shift between length and force. The cell length is im- 
posed as 

£=(£} + A e sin(ut), 

where {£) is the average length over an oscillation period, and 
Ae is the amplitude of the superimposed length-oscillations. In 
the linear regime, force is then given by 

F — (F) + A' F sin(urf) + A' F cos(urf). 

We now introduce the initial contact area An — tt(D/2) 2 , 
where D is the apparent cell diameter. The amplitudes are re- 
lated by 
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where 0' and 0" are differential viscoelastic stretching mod- 
uli. The storage modulus 0' reflects the purely elastic (non- 
dissipative) part of the cell reaction, and the loss modulus 0" 
the viscous (dissipative) contribution. In order to compare cells 
of different sizes and with biological gels, we use formal engi- 
neering stress units for the moduli. Accordingly, stress is taken 
as a = F/A . 

The differential stretch moduli 0' and 0" should not be con- 
fused with material parameters like the Young's modulus. The 
spatial distribution of force bearing structure inside the cell is 
unknown. Rather than introducing ad hoc hypotheses, such as 
assuming a uniform material, we simply treat the cell as a me- 
chanical black box. The unconventional symbol for the mod- 
uli intends to emphasize their experiment-specific nature. Fur- 
ther, these moduli characterize the response of the material to 
small perturbations around a situation which may be far away 
from the resting state. Indeed, we show below that it is suitable 
to study 0' and 0" as a function of the average force (F). A 
similar approach has been successfully used in stretching ex- 
periments on whole tissues, such as skin or muscle (uJ). Equiv- 
alent procedures have recently been applied to biopolymer gels 
under shear deformations, where the differential shear moduli 
are measured as a function of the average stress or strain 

Instead of the loss and storage moduli 0' and 0", it will be 
more convenient to regard the absolute modulus |0| and the 
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Figure 4: Lissajoux figures for different amplitudes. Top) Typical response at forces above crossover. Force F as a function of relative cell 
length £/{£) for strain amplitudes 3.5%, 6%, 13%, and 20%. The oscillation frequency is 0.2 Hz. Bottom) Typical response at compressive 
forces. Force Fas a function of relative cell length (./{£) for strain amplitudes 3%, 5%, 12%, and 20%. The oscillation frequency is 0.2 Hz. 



loss angle S, defined as 



|6| = 



arctan 



9" 2 
9' 



Fig -HI snows the dependence of the stiffness |0| as a func- 
tion of the strain amplitude A^/ (£). In general, at strain ampli- 
tudes in the range 0.02-0.06, the moduli do not change by more 
than 20%. Moreover, no significant distortion of the response 
is seen below relative deformations of 0.1, as shown in Fig.0] 
This holds irrespective of the frequency in the range 0.1-1 Hz. 
Thus, in subsequent experiments the amplitude is kept small, 
= 0.5 /im, which corresponds for all cell length values to 
0.02-0.06 strain amplitudes. More amplitude sweeps, includ- 
ing the amplitude dependence of the loss angle 5, can be found 
as supplementary material. 

For completeness, in Fig. we show preliminary results on 
the frequency dependence of the modulus 1 9 1 and the loss angle 
8, in the range 0.1-1 Hz. We are currently modifying the setup 
in order to explore a wider frequency range, as at present we are 
limited to just one frequency decade. Keeping this in mind, our 
results are consistent with the literature on cell mechanics. On 
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Figure 5: Amplitude sweep. Stiffness |0| as a function of the strain 
amplitude Ae/(£) for an arbitrary selection of cells. Each curve is a 
different experiment. All frequencies are 0.2 Hz. 
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Figure 6: Stress stiffening at constant length. The average cell length is kept constant throughout, (£) = 9^tm. Sinusoidal oscillations are 
superimposed with a strain amplitude Ai/{£) — 0.03. The frequency of the oscillations is cyclically changed from 0.1 to 1.0 Hz. a) The 
average force (F) is seen to increase in time, b) The modulus \Q\ increases in time for all frequencies, c) Stiffness \<d\ as a function of average 
force (F), for different frequencies. The line shows a power-law function y ~ x ' 7 . d) Loss angle S as a function of average force {F}, for 
different frequencies. 
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Figure 7: Frequency sweep (preliminary data). Modulus |0| and loss 
angle 8 as a function of frequency. Each curve is obtained using a 
different cell. 



Stress stiffening at constant cell length 

During the initial phase of force development after contact with 
the fibronectin-coated microplates, the cell sweeps force-space 
at a constant length. We superimpose sinusoidal oscillations 
to the constant average length, in order to probe the temporal 
evolution of the moduli 9' and 9". The frequency of the os- 
cillations is cyclically changed in the range 0.1 - 1.0 Hz. As 
shown in Fig. [6] as the average force increases with contractile 
activity of the cell, so does the modulus |9|. Figs.[6]c, d show 
the dependence of the response parameters |9| and 6 on the 
average force (F) for different frequencies. This is an exam- 
ple of active stress stiffening, since it takes place at an average 
constant length. 



the average, |9| increases weakly with the frequency, whereas 
5 is approximately constant. The modulus |9| increases with 
the frequency roughly as a power-law with exponents in the 
range 0.1-0.3. Such a frequency dependence, the signature of 
a flat, broad continuous spectrum of relaxation times, has been 
known in soft tissues for a long time I lb and has recently been 
observed at the micrometer scale 

Hi. 

Accordingly, recent 

experiments at the whole-cell scale performed with a setup sim- 
ilar to ours have revealed power-law creep functions Il2ri) . 



A master-relation characterizes stress stiffening 

Not all fibroblasts generate such high forces as in the experi- 
ment shown in Fig. [6] In order to span a larger range of aver- 
age force (F) and length (£), we step-stretch the cell and then 
keep the average length constant, superimposing oscillations to 
record the moduli 9' and 9". The procedure is periodically 
repeated, as shown in Fig. [8] As a reaction to a sudden change 
in length a force relaxation always occurs, usually followed by 
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Figure 8: The force response as result of imposed length changes. We 
step-strain the cell by about 50% at a rate of 1 .5 fim/s, and apply length 
oscillations at an amplitude of 0.5 fim and a frequency of 0.2 Hz. The 
experiment was performed using fibronectin mediated adhesion. 



active contraction. As the average force (F) evolves at a fixed 
length (£), the viscoelastic moduli are continuously recorded. 
We also perform step experiments controlling the average force 
(F). In this way, it is possible to span large areas in the (£) - 
(F) diagram. 

We see that both viscoelastic moduli depend only on the av- 
erage force, and are independent of the average length. The 
dependence of the loss angle 6 = arctan(8"/0') on the av- 
erage force (F) is erratic and weak, at most decreasing about 
20% in the whole force range. As a function of the individ- 
ual cell, it is within the range 10-30°. The absolute modulus 
|0| = (0' 2 + 0" 2 ) 1 / 2 remains constant at low forces, in a 
1-30 kPa range depending on the individual cell. Above a cell- 
dependent crossover force, we observe stress stiffening: |0| 
increases as a function of the average force (F). This depen- 
dence of |0| on the average force can be well approximated by 
a power-law, as shown in Fig.|5|f inset). More than one stress 
decade above crossover, most cells deform significantly and be- 
gin to detach or yield. 

A collapse of all data to a single master-relation can be 
achieved by introducing cell-dependent scaling factors, the zero 
force stiffness Oo and the crossover stress <Jc- On the average, 

( 9 for (er) < ere 




The exponent 7 is independent of the scaling factors. At 
26°C, 0.2 Hz, and 5% deformation amplitude, it is approxi- 
mately 1, as shown by the collapsed data in Fig.|9] The scaling 
factors are roughly related by 0o oc cr^.' 3 . Thus, an approxi- 
mate collapse can be reached with a single parameter. 
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Figure 9: The inset shows the elastic modulus |0| as a function of 
average stress for 13 cells, measured using length steps plus oscillation 
experiments such as in Fig.H0 The 

main plot shows the data scaled 
using 2 factors, which gives an exponent X ~ 1.0. All experiments 
are performed at 26° C and using fibronectin mediated adhesion. 



This master-relation is consistently found in all cells strong 
enough to reach average stresses above ~ 0.1 Qq. This repro- 
ducibility shows that oscillatory measurements are effective in 
probing cell mechanical properties independently of an under- 
lying slow active behavior. 

Ramp experiments with superimposed oscilla- 
tions 

The master-relation holds at a constant cell length while an 
underlying active contraction occurs. In this case, the rate of 
change of the average stress is set by the cell. As a means of 
assessing the validity of the master-relation while the average 
cell length changes, we increase (£) at a constant rate in the 
range 0.1-2 fim/s, and simultaneously superimpose small os- 
cillations at an amplitude ~ 0.5 /im and a frequency of 1 
Hz, to measure the dynamic moduli 0' and 0". As shown in 
Fig.^|fl, by stretching the cell, a change in average stress is in- 
duced. The average stress (er) depends roughly linearly on the 
average length (£) throughout a ramp. Remarkably, stress stiff- 
ening of the dynamic moduli is simultaneously observed. The 
master-relation between |0| and (er) is seen to remain valid at 
low deformation rates. These experiments show that the par- 
ticular way of sweeping force space is not relevant, since the 
|0|((cr)) function is qualitatively the same to that found in ac- 
tive contraction experiments. Only at rates higher than a cell- 
dependent value in the vicinity of 100 nm/s, |0|((cr}) falls be- 
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Figure 10: a) Ramp experiment with glutaraldehyde coating. Aver- 
age stress (top), loss angle (middle) and average length (bottom) as a 
function of time. Oscillations at 1 Hz, 5% amplitude are superposed 
throughout, also during the ramps (shaded areas). The phase differ- 
ence 8 increases with the deformation rate. 

b) Stiffness as a function of average stress. The relationship between 
|0 1 and (a) depends on the deformation rate. The curves correspond 
to the shaded areas 1, 2, 3 in Fig. llOl a. 

low the master-relation, and the cell becomes more fluid as evi- 
denced by an increase in the loss angle S (see Fig. [TO] a and b). 
In general, the master-relation cannot be deduced from the re- 
lationship between the average values (a) and (£), since these 
are linearly related in a ramp experiment. This absence of stiff- 
ening agrees well with results from ramp experiments on whole 
tissues (lullfjl) and fibroblast populated collagen gels (1151) . 

Cytoskeleton perturbation using drugs 

In order to distinguish between different cytoskeletal sub- 
components, experiments are performed in presence of drugs, 
which disrupt specific filament types. Only cells able to 
sustain tensions above ac are subsequently treated with 
drugs. We observe sudden detachment of the cells from 
glutaraldehyde-coated microplates with the actin depolymer- 



ization inducer Latrunculin-A at 1 /ig/ml iEtI) . The same ef- 
fect is seen with the myosin heavy chain ATPase inhibitor 
2,3-butanedionemonoxime at 2 mg/ml ilm . as well as with 
the myosin light chain kinase inhibitor ML-7 at 100 fjM J29I) . 
Since transmembrane proteins are covalently bound to the glu- 
taraldehyde coating, they must rip off from the membrane dur- 
ing drug-induced cell detachment. Thus, the cell membrane 
alone is not able to hold transmembrane proteins under signif- 
icant tension. An internal structure must bear the load under 
normal conditions. After disruption of either actin or myosin 
activity, this structure cannot sustain strong forces anymore. At 
4-fold lower concentrations, Latrunculin-A and ML-7 reduce 
|0| up to a factor of 5, while (a) goes to zero. No significant 
effect is seen with the microtubule-disrupting drug Nocodazol 
at 10 /ig/ml yOj). Taken together, these results show that the 
actomyosin system bears the tension, without any other signif- 
icant force-bearing structure in parallel. 

Stress stiffening with glutaraldehyde coatings 

In order to reduce active responses, we use glutaraldehyde- 
aminosilane coated walls, where accessible membrane pro- 
teins are covalently and non-specifically bound through imine- 
groups. Further, the serum concentration is reduced to 2%. In 
this way, the biochemical conditions are changed, but the ex- 
perimental cell geometry remains the same. Active responses 
are indeed found to be weaker in these conditions. Compared 
to fibronectin mediated adhesion, the stiffness values are scaled 
down by about a factor of three. However, the master-relation 
between |0| and (a) can still be observed. 

DISCUSSION 

The master-relation is a general feature, which reproducibly 
occurs in our experiments. It holds during active contraction 
and adaptive responses, such as those seen in Figs.|5J |U and 
[8] as well as in ramp experiments with superimposed oscilla- 
tions, as in Fig.^|fe. It is observed in different biochemical 
environments, which change the intensity of active responses. 
From this generality of the master-relation, we conclude that 
it reflects a fundamental property of the cell force-bearing ele- 
ments. 

Stress stiffening in actin networks 

Our cytoskeleton perturbation experiments point towards the 
actin network as the main component responsible for stress 
stiffening. Actin filaments, with a persistence length of 15 /im 
QEl), behave as semiflexible filaments in vivo, where typical 
filament lengths are about 1 /im. It has recently become ap- 
parent that crosslinked as well as entangled networks of semi- 
flexible filaments show characteristic nonlinear mechanical be- 
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havior. In vitro prestressed biopolymer gels under shear defor- 
mation, including crosslinked actin networks, generally show 
a transition from a linear regime to power-law strain stiffen- 
ing (0 0, [32I [33I [34l> . The connection to our results has to be 
established with care. In vitro gels are passive minimal sys- 
tems, generally studied under shear deformations. Instead, the 
living cell is a complex entity where a multitude of crosslink- 
ing proteins are available, many of them of a dynamic nature. 
Cytoskeletal restructuring might take place within an oscilla- 
tion period. Keeping this in mind, we still find the qualitative 
similarity between the stiffness-force relations very suggestive. 
In scruin-crosslinked actin networks (|7|), the ratio between the 
crossover stress and the zero force modulus is of the order of 
10 -1 , and the stress-stiffening exponent is in the range 1-1.5, 
in remarkable agreement with our results. In experiments on 
filamin-crosslinked actin gels (0), where the storage modulus 
G is actually studied as a function of strain, the crossover strain 
is also of this order. 

Intrinsic stress stiffening 

The master-relation connects the viscoelastic moduli to the av- 
erage force independently of cell length. Rather than strain 
stiffening, fibroblasts exhibit intrinsic stress stiffening. To rec- 
oncile this to our interpretation of the master-relation, we pos- 
tulate that the internal strain of the gel, i.e., the stress, uncou- 
ples from the cell length £ as a result of active contraction, via 
rearrangements of the connection points of the network. 

As shown by ramp experiments with superimposed oscilla- 
tions, the master-relation holds even during an externally im- 
posed length increase, if the extension rate is sufficiently slow. 
The magnitude of the deformations which can then be imposed 
(more than 50%), without the cell yielding, suggests that ad- 
dition of material takes place. As the loss angle 6 does not 
change under such conditions, we speculate that the cytoskele- 
ton grows without major structural changes. The cytoskeleton 
is certainly capable of such restructuring, as exemplified b y th e 
growth under compression of the actin tail of Listeria 
towed growth in neurites l3m or stress fibers losing material 
while contracting TO . The maximal growth speed from our 
ramp experiments would be around 6 /im/min, in good agree- 
ment with the growth of the actin tails of the bacteria Liste- 
ria 10. At higher deformation rates the cell starts to fluidize, 
which can be interpreted by an increased fraction of ruptured 
bonds. 

Master-relation in the context of cell mechanics 

Several features of the master-relation can be observed in ex- 
periments performed at very different length scales, corroborat- 
ing our interpretation as a general feature of the force-bearing 
elements in the cell. This shows the master-relation to be gen- 



eral, independent of the type of mechanical experiment or cell 
geometry. Uniaxial stretc hing exp eriments performed on skin, 
myosin fibers, tendons qjj, iKj Il7l) and on fibroblast-populated 
collagen gels ( 15) show proportionality between dynamic stiff- 
ness and force. These results can be seen as a particular case of 
power-law stress stiffening, for an exponent 7=1. 

Our results are also similar to stress stiffening in skeletal 
and smooth muscle, where stiffness is proportional to force. 
The generally accepted explanation is that both stiffness and 
force are a function of the variable number of actomyosin cross- 
bridges 13% . Although such an explanation is attractive, it does 
not seem to apply to our case. We have measured forces up 
to 1 /iN, in agreement with total forces exerted by spreading 
fibroblasts of ~ 2 /iN (0 H3). This corresponds to - 10 6 
myosins working in parallel, very close to the total amount re- 
ported in fibroblasts Il4lh . Taking crossbridge stiffness as 0.6 
pN/nm (13), an arrangement of 10 6 myosins in parallel would 
be a factor of 100 stiffer than the maximum |0| we have mea- 
sured in fibroblasts. 

At the subcellular scale, microrheology experiments per- 
formed on adhering cells show proportionality between the 
shear storage modulus G' of the actin cortex and the force ap- 
plied by the cell on the substrate Here, the force was 
increased by stimulating cell contractility with histamine, or 
decreased by the relaxing agonist isoproterenol. The result 
compares well to our observation of stress stiffening at con- 
stant length, though no crossover to a linear regime is reported 
here. Recently, simultaneous increase of both the storage and 
loss shear moduli G' and G" of the actin cortex after stretch- 
ing epithelial cells has been observed Accordingly, we 
see stress stiffening of both longitudinal moduli 0' and 0" 
when stretching the cell in a ramp experiment with superim- 
posed oscillations. Thus, the fact that the master-relation holds 
regardless of the way force-space is explored agrees with these 
microrheological studies taken together. 

A similar example of length scale invariance is given by the 
agreement between microrheological studies in the frequency 
domai n Jl4l E5I) and whole-cell experiments in the time do- 
main (26). Fabry et al studying the shear moduli G' 
and G" of the actin cortex as a function of frequency, could 
collapse data from cells in different conditions onto master 
curves. The response in the frequency domain as revealed by 
these microrheological procedures was subsequently found to 
agree with the power-law creep function of single cells H^ti . 
Hence, scale invariance and self similarity may well be gen- 
eral properties of cell mechanics in the force, as well as in the 
time-frequency domain. 
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Possible stiffening mechanisms 

Aiming at describing crosslinked biopolymer networks, theo- 
retical frameworks have been developed recently. Networks are 
modeled as a 2-dimensional random arrangement of filaments, 
characterized by a bending and a stretching modulus. Tem- 
perature effects have been addressed by including an entropic 
stretching modulus At a non-zero temperature, thermal 

energy is stored in bending fluctuations, so that the distance 
between crosslinks is smaller than the actual contour length of 
the filament. Separating these crosslinks by stretching reduces 
the amplitude of these shape fluctuations, which gives rise to 
an increased effective stretching modulus. In general, the me- 
chanical response of 2-dimensional random networks depends 
strongly on the length scales involved, the filament length and 
the mesh size £ (the average distance between filaments). For 
long filament lengths or small mesh sizes, the deformation of 
the network is affine, i.e., the macroscopic mechanical response 
is essentially given by that of single filaments (|8|). The non- 
linear elasticity of crosslinked biopolymer networks has been 
explained in terms of random networks (0 0), as an affine 
stretching-dominated regime where the macroscopic response 
is given by single-filament entropic stretching. This response is 
highly nonlinear, since thermal fluctuations are suppressed as 
the filament ends are separated close to full contour length. At 
high forces, the stiffness-force relation becomes a power-law 
with exponent 3/2 |7J). 

An alternative explanation for stiffening in crosslinked ran- 
dom networks has been recently presented by Onck et al 
Here, stiffening arises as a transition from bending to stretch- 
ing. Thermal effects are mostly irrelevant; increasing the tem- 
perature from to 300 K increases the crossover strain, but does 
not affect the stiffening regime. Also the transition to stiffening 
arises in spite of a certain degree of non-affinity. 

Below we point out some remarkable features of single fila- 
ment bending response which may be relevant for crosslinked 
filament networks. We discuss the bending response of inex- 
tensible filaments, characterized by a length L and a bending 
modulus k. We show that bending itself shows a crossover to 
power-law stress stiffening with an exponent 1.75, as well as 
a crossover strain of the order of 10 _1 . This can be seen in 
Fig At large forces, the stiffness-force relation becomes a 
power-law with exponent 3/2, similar to the entropic stretching 
response. The details of the calculation are left to the appendix. 
It is noteworthy that many features of the master-relation are 
similar to the mechanical response of single filaments. The ex- 
perimental fact that a single parameter is sufficient to obtain 
the master relation may be captured by both the force and stiff- 
ness scales varying as k/L 2 . Interestingly, the magnitude of 
the experimental crossover stress, when expressed as a strain 
ac/Qo, is of the order of 10 _1 , in good agreement with the 
theoretical crossover strain. As Fig.^Jfr shows, the first force 
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Figure 11: a) Single-filament nonlinear elasticity: filaments are re- 
garded as inextensible, characterized by a contour length L and a 
bending modulus k. One end is clamped, the other one free. The 
force F bends the filament by an amount x. The dimensionless force 
is defined as / = FL 2 /k. 

b) Differential stiffness df/d(x/L) as a function of the dimensionless 
force /. 

c) Euler-Beraoulli elastica: filament shapes for different forces. 



decade above crossover of the theoretical stiffness-force rela- 
tion is an approximate power-law with exponent 1.7, within the 
experimentally observed range. Further, the magnitude of the 
force scale is the right one if one assumes a realistic cytoskele- 
ton mesh size of 100 nm ( SHE! and an actin bending mod- 
ulus of k = 60 nN(nm) 2 Q. This corresponds to the open- 



cell foam model proposed by Satcher and Dewey as a general 
model for the cytoskeleton (I44I) . which gives a zero-force stiff- 
ness 0o ~ 10 kPa, in goo d agreement with our measurements 
and with literature Jl7l t44h . 

These observations suggest that the bending response of fil- 
aments may play a role in the nonlinear mechanical properties 
of crosslinked gels. Therefore, it is worthwhile considering this 
effect in further modeling of such systems. With the present 
set of experiments, we cannot distinguish between the differ- 
ent mechanisms proposed for stiffening. Further experimental 
work and theoretical modeling will be required to elucidate the 
relative contributions of these effects for physiologically rele- 
vant values of mesh sizes and filament lengths. 



Summary, conclusion and outlook 

By performing single fibroblast mechanical measurements we 
have revealed a master-relation relating stress and dynamical 
cell stiffness. This relation is obtained by a simple scaling of 
data from different cells. For the measurement we use the fact 
that cells are active: as a reaction to a sudden perturbation the 
cell sweeps a range of mechanical stresses. We probe the cell 
elastic properties during this active response. In spite of this 
underlying complex behavior which may involve multiple bio- 
chemical pathways, the master relation is surprisingly simple 
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and reproducible. If the average force is externally changed, 
by slowly stretching the cell in a ramp experiment while su- 
perimposing small-amplitude oscillations, the master-relation 
between (F) and the stiffness |0| is seen again. Thus, a dis- 
tinction between active and passive stress is artificial - regard- 
less of deformation history, the response to small perturbations 
is always given by the average force. In view of this general- 
ity, and the remarkable agreement with the mechanic behavior 
of biopolymer networks, we interpret the master-relation as re- 
vealing the nonlinear response of the actin network. 

These experiments provide a framework where simplified 
biomimetic systems can be connected to living cell mechan- 
ics. Here, the active nature of the cytoskeleton is naturally in- 
tegrated into the phenomenology. Along these lines, we are 
currently extending our frequency and amplitude range, to un- 
dertake a more detailed analysis of the master-relation. 

As a final remark, fibroblasts reach the stress-stiffening 
regime naturally by active contraction and we expect the 
crossover force to be within the physiological force-range. 
Thus stress stiffening may play a role in vivo. For example, 
cells are known to sense the extracellular elasticity and exert 
forces along lines of maximum effective stiffness, a process 
which has been modeled within linearized elasticity yj). In 
the future, power-law stress stiffening could be included in the 
modeling of cell organization by mechanosensing. 
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F. Jtilicher and R. Everaers for helpful discussions, and B. Wehrle-Haller for his 
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care of the cell culture is gratefully acknowledged. We also thank Y. Marcy 
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reviewers for constructive criticism and helpful suggestions, and to the editor 
for his interest in the manuscript. This work has been funded by the Universitat 
Bayreuth and the Institut Curie (Paris). Some of the results of this work were 
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and in the DPG Fruhjahrtagung, Regensburg, Germany, March 2004. 

Appendix 

Bending response of an inextensible filament 

In what follows we show that a transition to power law stiff- 
ening is already present in the nonlinearity of bending defor- 
mations. We analyze the mechanical response of inextensible 
filaments characterized by a bending modulus n and a contour 
length L. We assume that the relationship between bending 
moment and curvature stays linear, so that any nonlinearity is 
of geometrical nature. The filament shape then corresponds 
to Euler-Bernoulli elastica d46l) . Our boundary conditions are 
one filament end free and the other one clamped perpendicular 
to the force direction, as illustrated in Fig .^J a. These con- 



ditions are the minimum requirements for bending. At small 
strains x/L, the mechanical equilibrium equation can be lin- 
earized and the zero-force stiffness analytically calculated as 
dF/dx = 3k/L 3 . 

Numerically solving the nonlinear equilibrium equation, we 
find a relatively sharp crossover to a strain-stiffening regime at 
strains x/L > 0.3, as shown in Fig.[^fr. At large strains, the 
stiffness-force relation asymptotically becomes a power-law 
with an exponent 1.5. At high forces, the filament shape tends 
to a straight line along the force direction, with a sharp kink 
at the clamped end (see Fig.[^c). In this limit, the effect of 
increasing F can be shown to be a scaling of the filament shape 
at the clamped end; the endpoint coordinates and curvature ra- 
dius remain linearly related as they change as a function of 
F. Thus, we can introduce a force-dependent scaling factor e, 
which tends to zero as the force is increased to infinity. In terms 
of this factor, the full contour length L and the projection length 
x (see Fig.[^a) can be written as L — x oc e, and the torque 
balance equation as Fe oc 1/e. Therefore, F oc 1/(L — x) 2 , 
which gives power-law stiffening, dF/dx oc F 3 / 2 . 

In general, an entropic contribution to the longitudinal fil- 
ament stretching modulus is expected due to thermal fluctua- 
tions. At low forces, this leads to a thermal longitudinal modu- 
lus k T = 90K 2 /(fc B TL 4 ) 0), where L is the contour length, 
and ksT the thermal energy. Taking a mesh size of 100 nm, this 
thermal longitudinal modulus fc^ is about a factor of 10 4 larger 
than the low-force mechanical bending stiffness, k± = 3k/ L 3 . 
Since k±/kx oc L, the mechanical response of a regular net- 
work such as the open-cell foam model (44) remains dominated 
by bending stiffness for physiological mesh sizes. Further, for 
such mesh sizes the mechanical bending energy already ex- 
ceeds the thermal energy at strains above ~ 4%, well within 
the linear regime. Thus, thermal effects are not relevant for the 
bending response itself. 
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